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Abstract 

We develop a systematic multi- local expansion of the Polchinski- Wilson exact renor- 
malization group (ERG) equation. Integrating out explicitly the non local interactions, 
we reduce the ERG equation obeyed by the full interaction functional to a flow equa- 
tion for a function, its local part. This is done perturbatively around fixed points, but 
exactly to any given order in the local part. It is thus controlled, at variance with 
projection methods, e.g. derivative expansions or local potential approximations. Our 
method is well-suited to problems such as the pinning of disordered elastic systems, 
previously described via functional renormalization group (FRG) approach based on 
a hard cutoff scheme. Since it involves arbitrary cutoff functions, we explicitly verify 
universality to 0(e = 4 — D), both of the T = FRG equation and of correlations. 
Extension to finite temperature T yields the finite size (L) susceptibility fluctuations 
characterizing mesoscopic behaviour (Ax) 2 ~ L 9 /T, where 9 is the energy exponent. 
Finally, we obtain the universal scaling function to C(e 1//3 ) which describes the ground 
state of a domain wall in a random field confined by a field gradient, compare with 
exact results and variational method. Explicit two loop exact RG equations are derived 
and the application to the FRG problem is sketched. 
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1 Introduction 



The idea of writing an exact equation for the scale dependence of the full action functional 
already appears in the review of Wilson and Kogut 0. Since it is an equation for a full 
functional of the fields, its detailed analysis is hindered by technical complications. The 
much simpler Wilson momentum shell integration method is commonly used for one loop 
calculations. Since it does not follow the full functional, subsequent efforts were made to 
embed it into a better controlled sharp cutoff exact RG [Q, Qj. For practical perturbative 
calculations beyond one loop, field theoretical renormalization methods are more often used 
since they have proved vastly more efficient. However, the exact RG equations offer the hope 
to develop ab initio calculation relying on no assumption, possibly non perturbative, from 
any bare model. In principle it should be useful to obtain precise results when applied to bare 
theories for which we have little insight on possible underlying field theoretical description. 

In the work of Polchinski jlj, the exact RG equation was put on a more precise and 
aesthetic framework, and used to prove the renormalizability of the 4 theory in four dimen- 
sions. The exact renormalization group equations indeed provide formal results or general 
proofs about symmetries ||. For practical calculations however, one needs to truncate in 
some way these highly complicated functional equations. To do so, different procedures 
have been proposed || 0], and have been mainly applied to the study of non-perturbative 
problems ||. For example the exponents of the 0(n) model in three dimensions were es- 
timated 0] using a choice of truncation. One commonly used projection method is the so 



called local potential approximation ||10|| , obtained by a constant background field method 
neglecting the momentum dependence. Further extensions include additional projections on 
higher gradients of the field |TT[| . Although very interesting, these projection methods are 
often uncontrolled. More accurate results are expected if more couplings are kept, which is 
possible with heavy numerical integrations of flow equations. In this respect, exact RG as a 
tool is now used both in particle and condensed matter physics. For instance, outstanding 
problems in strongly correlated electrons such as the Hubbard model in D = 2, have been 



recently studied by numerically integrating the flow of a large number of vertices |T2|, using 
a fermionic version of the Polchinski equation [|J. 

By contrast, comparatively a few works use exact RG method to develop perturbative 
calculations. One example is the computation of the beta function of 4 in 4— e dimensions to 



one loop|13[, where universality is made particularly explicit through the use of an arbitrary 
cutoff function. Although obviously more powerful methods are available in that case, there 
are some problems in condensed matter physics which appear within reach of perturbative 
calculations but for which no coherent field theoretical formulation is available at present. 
This is the case for the pinning of an elastic system in a random potential, for which a 
momentum shell RG method has been developed [|14|]. In this problem, an infinite number 
of coupling constants becomes relevant for D < D c = 4 and one must write a RG equation 
for a full function R(u) (the second cumulant of the disorder), hence the name "functional" 
renormalization group (FRG). As such it differs from standard field theoretical RG. Thus, to 
understand better this problem, i.e to show explicitly universality to one loop and beyond, 
there is a need for a perturbatively controlled exact RG method, able to admit a full function, 
the local part R(u) as a small parameter. Indeed the field theoretical formulation is frought 
with difficulties, in particular because the function R{u) develops non analytic behaviour at 
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finite scale. These issues are discussed in a related work [[OJ. 

In this paper we develop a novel method to solve the Polchinski exact renormalization 
group equation and use it for explicit calculations. Writing the action as a sum of multilocal 
interactions, we note that the Polchinski equation naturally reduces to a hierarchy of equa- 
tions obeyed by simple functions. This hierarchy can be solved in an expansion in powers of 
the local part. Indeed, we find that exact integration of the multilocal parts yields a single 
RG equation for the local part. The method is thus controlled around fixed points where 
the local part is proportional to a well defined small parameter (e.g. e = 4 — D). It does 
not require any arbitrary projection procedure or neglect of operators, as is usually done 
in derivative expansions or local potential approximations. In addition, we obtain explicit 
formulas for any correlation function which allow for practical calculations. Since this is 
done for arbitrary cutoff functions, it allows explicit check of universality order by order in 
the expansion. 

The aim of this paper is twofold. On the one hand we present the general method to 
all orders, valid for a large class of theories. We derive the explicit form of the exact RG 
equation for the local interaction up to third order. On the other hand we apply this method 
to several problems, first as a check, to the 0(n) model, and second to the FRG for disordered 
elastic systems. Explicit calculations and applications in this paper are restricted mainly to 
one loop. Although briefly mentionned here, applications to two loops will be detailed in a 
companion paper |16] . 



Two variants of the method are presented. The most direct one consists in a straight 
expansion of the action in multilocal terms. The second one consists instead of first absorbing 
tadpoles into the interaction (so-called Wick ordering), then expanding. Being inequivalent, 
they provide independent checks of the universal results. The first method yields more 
complicated equations but, can be better suited to some problems, such as the T > FRG. 
Note that although Wick ordered versions of the Polchinski equation have been studied 
before, the multilocal expansion performed here is to our knowledge novel. 

As mentionned above, the method is indeed well-suited to the FRG for disordered elastic 
systems of internal dimension D since there the full local part is controlled by e = 4 — D. 
It allows us to show that the one loop FRG equation, as well as correlation functions, 
are independent of the cutoff function. In addition we obtain higher cumulants of the 
renormalized disorder, which as the second cumulant, are non analytic functions. This is 



necessary to escape the so-called dimensional reduction i.e the property of the present 
theory by which all perturbative calculations at T = are identical to the same calculation 
in a trivial gaussian theory |T7| (see Appendix |IJ) . This nonanalytic behaviour is rounded at 
finite temperature T and we obtain the scaling form of the rounding region. This allows us 
to compute, for the first time using the FRG method, the susceptibility fluctuations which 
characterize the glassy behaviour of finite size systems. Finally, we obtain the universal 
0(e 1 ^ 3 ) correlation function which describes the ground state of a domain wall in a random 
field confined by a field gradient, compare with exact results and variational method. 

The method presented in this paper also allows to investigate the theory of disordered 
elastic systems beyond lowest order in e (one loop). A recent two loop calculation was 
presented in [j44] |. However since it was performed at T = 0, and for analytic R(u) it 
fails beyond a finite length (Larkin length) and cannot describe universal properties. The 



application of exact FRG to next order is described in [16]. We sketch here however some 
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preliminary results. 

The paper is organized as follows. In Section ^J] we present in a pedagogical way the 
conventional exact RG method. Appendix |A| and Appendix [B] provide complements, respec- 
tively about general invariance properties of the correlations and about examples of solvable 
cases of the Polchinski equation. In 2^2, the multilocal expansion in the local part U is 
introduced up to bilocal terms. The ensuing RG equations to order U 2 are given in |273|. 
The multilocal expansion to arbitrary order and the RG equation to order U 3 is given in 
Appendix y. The multilocal expansion of the Wick ordered functional up to second order, 



and the resulting one loop RG equation is presented in |3J]. The general multilocal expansion 



(Appendix P.l|) and the resulting form to third order is given in Appendix G2 . The explicit 
two loop RG equation is obtained in |G.3| . Application to the 0{n) model to one loop is 
presented in jO[ We then turn to applications to disordered elastic systems in Section f| 
(one loop) and Section ^| (two loops). First we recall and generalize in Appendix |B| the 
dimensional reduction phenomenon. Then the T = FRG equations are established in Sec- 



tion ^?T| and finite temperature extension are given in Section [4.2| . Finally, the calculations 
of the scaling function in the random field Ising model is performed in [4.3| . We sketch in 
Appendix [E] the variational calculation to be compared with the FRG results of O and 
sketch some preliminary steps of a two loop FRG in Appendix [F]. 



2 Method 

2.1 Exact RG procedure 

Consider a system whose state is described by a bosonic field 0f = (j>i(x), where x denotes 
position in space, and i is a general label denoting e.g. fields indices, spin, replica indices, 
additional coordinate (e.g. time) etc. (or more generally any quantity which will not undergo 
the coarse-graining). The system, in the presence of external sources Jf , is described by the 
partition function: 

Z{J) = [ e J ^~ sw (1) 

J 

obtained by the integration over the field 0, where the action <S(0) is a functional of the 
field 0, and J : denotes here and in the following the full scalar product (e.g. f x J2i Jf 4>f > 
with J x = / d D x). In a problem of equilibrium statistical mechanics, S(4>) = Ti,((/))/T where 
T~i{4>) is the hamiltonian and T the thermodynamic temperature, the free energy being 
F = — TlnZ(O). Averages of any observable A((f>) (i.e. functional of 0) are defined by 

(A<f>))s = j-^im ( 2 ) 

The usual way to compute correlation functions and averages is to perform a perturbation 
expansion, writing the action as a sum of a quadratic part and a non-linear part V(0) 

5(0) = ^0 : G- 1 : + V(0) (3) 
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where V(4>) a functional of 



and G^f 



In the following we denote the Gaussian average of any observable by 



^-^ij Ixy 1 

[A(4>)]g = (-4(</ ) ))s G with respect to the quadratic theory Sq 
the generating function of all correlation functions 



is a symmetric invertible matrix, : G 



: G 



-i 



W(J) = In 



G 



We introduce 



(4) 



Note that it differs from the usual definition by a J-independent quantity lnZ(J) = W(J) + 
|TrlnG. The ultra-violet cutoff, present in physical models, is necessary to yield finite 
results in the perturbative calculation with respect to V. A broad class of soft cutoffs can 
be implemented on the Gaussian part, giving a vanishing weight to fast fields. For example 
a scalar massive theory, rotationnally invariant, is regularized in the UV by the following 
general cutoff function 



G(q) 



• q ' 

■2A 2 - 



q 2 + m? 



(5) 



where c(0) = 1 and c(s) decreases rapidly to zero for s > 1 as in Figure [I]. 



c(x) 




A' A 

Figure 1: Cutoff function. 



Ao 



The exact RG method |], |J consists in varying the cutoff A and writing an equation for 
the function V(0) so as to conserve exactly the averages of all observables involving only 
"slow" modes of the field. More precisely, the average of an observable A(<f)) depending only 
on modes q < A' of the field can be computed within any of the theories linked by the 
equation presented by Polchinski corresponding to a cutoff A > A'. As in any RG procedure, 
the strategy will be to compute averages of slow observables using the coarse-grained theory 
of cutoff A' < A . 

To this aim, a set of actions Si 



5 I (0) = -0:Gf 1 :0 + V,(0) 



(6) 



is introduced, where the Gaussian part is an arbitrary function Gi of I (e.g. corresponding 
to a cutoff A; = A e - '). The initial propagator, corresponding to a cutoff A , is denoted 
Gi=o = G and the bare interaction V;=o = V. The correlation functions in Si derive from 
Wi(J) = In \ e J--<l>-ViW] = \ n Zi(J) - ^TrlnG,. 

L J Gi 

For any given A', one defines a (A'-)slow observable to be a functional of <fi depending only 
on the 4> q = (j)(q) with q < A' = Aoe~ l . We want to choose the /-dependent non-linear part 
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Vi(4>) so that the averages of slow observables remain unchanged. Through differentiation, 
it is equivalent to ensure that 

diWi(J) independent of J (7) 
for any source J with J q = for q > A'. Using the general identity 

dilMfiU = ^Tr (d^ : l-^Mtih) + [diA^)] Gl (8) 

valid for Gaussian averages of any /-dependent observable Ai(<p) and applying it to Ai{(p) = 
e J*h-vM) ^ one nnc i s 



\{J- ^ W)) : d l G l : (J - A W )) - d l V l {ct>) - l -Tvd l G l : ^V(<f>)) e**-*® 



G, 



where here and in (f|), TrA : B = J xy 'Bj* . Hence, if Vi(4>) satisfies the Polchinski 
functional equation 

mW = Ug, : *L W) ) + i ' W ) : m : j-V M ) (9) 



then the above conservation condition ([FD is satisfied. We have used explicitely the condition 

J : did = (10) 

which imposes that the cutoff function verifies diGf = for q < A' and I > I'. Hence, for the 
example fl5|) one has to choose [19[ cutoff functions c(s) such that c(s) = 1 for < s < so 
with some (arbitrary) s . 

The above framework is in fact too restrictive. We can easily lift the restriction on slow 
modes (and on the form of the cutoff function c(s)). The applications of Polchinski equation 
can be generalized to the computation of any observable (not restricted to be "slow" ) . As 
shown in Appendix |A.1| one can indeed express W(J) in terms of any of the /-dependent 
actions Si(<f>): 

W{J) = ^J:(G-G:Gr 1 :G):J + W t {J : G : Gj 1 ) (11) 



In fact we show in Appendix [A.2| an even more general method which allows for arbitrary 
field rescalings. 

Differentiating W(J) once yields ((f)) s = G : GY 1 ((f))s l , and once again yields the two 
point connected correlation function 

(<M>)% = G + G:Gi 1 : (<#>^ - G t ) : G^ : G (12) 

and so on for higher correlations. When performing a perturbative calculation, the factors 
G : G^ 1 restore the original propagator for the external lines, whereas internal lines of the 
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graph involve Gi. Accordingly, with this procedure the function c(s) can be arbitrary (it is 
however convenient - see below - to use c'(0) = 0). 

Note that if G : J = Gi : J, one recovers (^), i.e. W(J) = Wi(J) for these slow J's as a 
special case of (|). In that case, for q < A', (|T2| ) reduces to (00)s = {<P4>)si as ^ should. 

To compute correlation functions, it is useful to express W ( J) in a perturbation expansion 
in powers of V;(0), which reads to lowest order 



W{J) = ~J:G:J-e 



J:G:G~ L :G:.J 



e J:G:G~ 



J Gi 



(13) 



The Polchinski equation @ can equivalently be written as a functional "diffusion" equa- 



tion 



-V,(0) 



or in its integrated form 



-Vj(0) 



-V (</>+</>') 



G-G 



(14) 



(15) 



where the average is over 0', which makes explicit the definition of Vi(<p) as a coarse-grained 
interaction, i.e. integrated over the "fast part" <$' of the field. In fact, the decomposition into 
slow and fast modes and the definition of coarse-grained observables relies on the property 



[A(<P)] G = [A{4> + 



J G-G, 



Gi 



of Gaussian averages (see Appendix |A.1[ ) . 



Although in general the Polchinski equation is far too complicated to be solved, in some 
simple cases one can find exact solutions e.g. Gaussian models, zero dimensional toy model. 
Most interestingly, there exist a large class of exact solutions which appear as superpositions 
of gaussians. In all these cases, one can explicitly verify an interesting property of the 
Polchinski equation to generate cusp singularities. This is further discussed in Appendix |TJ 
and in forthcoming publication [05 . 



2.2 Multi-local expansion 

The Polchinski equation, in addition to being elegant, is conceptually more satisfactory 
than other RG methods, e.g. Wilson's shell renormalization, because it is exact and better 
controlled since. Being valid for arbitrary cutoff procedures, it does not suffer from the 
problems associated with the sharp cutoff 0. However, this functional equation generates 
non-local operators, which until now, has limited its practical applications. This generation 
can be seen in terms of Feynman diagrams and compared to Wilson's shell renormalization, 
since G — Gi which contains a range of wave- vectors centered around Aoe _z , plays the role 
of the on-shell propagator. The term V" with a second derivative in (|j) represents tadpoles 
while the term V'V represents diagrams with only one contraction (one particle reducible). 
These last terms are non-local operators. For instance in 4 theory, it generates the operator 
4>(x) 3 dG x ~ y <f)(y) 3 which is bi-local since it corresponds to a graph where external momenta 
must be greater than A;. The way Polchinski's equation reproduces the loop diagrams 
(i.e. local terms) is that after integration over a slice dl, a bilocal interaction generated by 
the second term of ® is fed into tadpole diagrams. A fast momentum goes around the 
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corresponding loop, and slow external momenta are allowed. Thus one needs to integrate 
the flow and study the feedback of the generated non-local operators into local ones. 

We now present a method which allows to perform this program in a controlled way. The 
following expansion in the number of points (local, bi-local, etc..) 

V(0) = / U((f> x ) + / V(F, <f>y,x-y) + ... (16) 

J x J xy 

is valid a priori for any translationally invariant functional V(<f>) interactions. We discuss 
here only the first two terms, the general systematics being given in Appendix [Cj Here, U((f>) 
is a function of the vector fa and involves the value of the field at one point in space. The 
bi-local part is a function V(4>, if), z) of two vectors <p, if) and a space coordinate difference z. 
In order that the expansion be well-defined, one needs the bi-local interactions to have no 
projection on the local ones. A natural way to define such a projection, inspired from the 
conventional short-distance-expansion, is the exact equality 

/ F(r,<P y ,x-y) = f F(</?,F,y)+ f (P(0^ 2 , 0^/ 2 , y) - F(f, <fT, y)) (17) 

Jxy J xy J xy v ' 

= f(F 1 F)m+ [ {{l-P 1 )F){<jfA v ^-v) (18) 

J x J xy 

where we have introduced the projections 

(P x F)(0) = / Fifi, fry) (19) 
(P 1 F)(<f>,iP,z) = 6(z) I F((f),if),y) (20) 

on the subspaces of local and bi-local interactions respectively. Indeed, (Pi(l — Pi) F) ((f)) = 
J x ((l — Pi)F)((f),(f),x) = 0, i.e. (1 — Pi)F has zero local part and is thus properly bi-local. 
Interestingly, this definition implies that the function V((f),if),z) appearing in the proper bi- 
local operator of ([H]) also satisfy the stronger property J z V((f>,if>,z) = for any <f>,if). Note 
that with no loss of generality, V(<p, if), z) = V(ip, (p, —z). Here in addition, we will consider 
parity invariant theories (V((f>,ifj, z) = V((f),if>, —z) too). 

For theories where the initial interaction U is local and is formally treated as a "small" 
quantity U (e.g. the 4 theory in D = 4 — e where U ~ e), it is natural to consider that the 
bi-local term will be of higher order 0(U 2 ). In fact, this property results from the Polchinski 
equation since the term which creates bilocal interactions from local ones is 0(U 2 ) (the first 
part V" does not increase the degree of non- locality of V). This property that solutions 
of the Polchinski equation can be organized in powers of U depending on their locality 
holds to arbitrary orders (p-local operators are 0(U P )) as is discussed below and shown in 
Appendix y. Thus, to lowest non-trivial order 0(U 2 ), the flow equations involve local and 
bilocal parts. Their schematic structure is 

dU = U" + Pi(V" + U'U') + 0(UV) + 0(V 2 ) (21) 
dV = (1 -Pi)(V" + U'U') + 0{UV) + 0(V 2 ) (22) 

where we have written the subdominant terms which will be neglected in the following. 
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A simplification occurs if we choose, as done in this Section, (diGi)ij(q = 0) = 0. Indeed, 
the term P\VV vanishes since 



frGf^jdiUMdjUte) = (23) 
Let us write now in an explicit form: 

Wii!l>) = —dG^d&Uticj)) - J^dld'V^^x) (24) 

W(<f>,1>,x) = -did] (d(^V l (M,x)-6(x)Jd(^V l (M,y)^ (25) 

-\dG%* (d]d] + dtdf) Vfa^x) + IdGZMWdjUM) (26) 
where dG stands for diGi and d}A{<p,ip) (resp. dfA(<fi, ip)) for -^pA((p,ip) (resp. -£pA((f), ip)). 

2.3 Solution to the lowest order RG-equations 

To solve (21j22|), one switches to Fourier space (in the field): 

U K = J d<j)e- iK ^U{<j)) (27) 

V KPx = [ d<j)d^e- iK ^- iP ^V{<j),ifj,x) (28) 



where K.(f> = Yl,iKi$i- It turns out that the equation for V\ can be integrated explicitely as 
a retarded function of Iff. 

vr x = l{ F i KPx - 6 ^J y F i KPv ) ( 29 ) 

Jo 



since we have chosen V^ Px = 0. One can then reinject this result in (21) and obtain a closed 
RG equation for Ui(<p): 

diU^) = —dG^didiUi^)-]-! e^ K+p ^ f K.(dGf-dGr°)-P 
l l Jkp Jx 

t dl' K.dGl.Pe^ K - {G r°- G T ) ■K + IP.(G f =°-G r ).P+K.(G f -G tf ).P u K pP ^ 

Jo 

This is the exact renormalization equation for an arbitrary local interaction U(<f>) to 0(U 2 ). 
Note that the second order term is retarded in I, since as discussed above, local terms are 
generated only after integration. 

More generally, this procedure can be carried out to any order in U using the hierarchical 
structure of the flow equations for p- local interactions (see Appendix It is found that 
the general structure for the flow of the local part is 

dflitt) = P[Ui><i{<!>')]{<t>) (32) 

= Ejf , , ^.^[^-^nM-'-U^l^ (33) 

n>l Jh<h<-<l n <l 0<pi 0(p n 
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and in ( |3"TD we achieved the calculation of the (3 function to second order in U. Once the 
solution of (|32|) is known up to 0(U P ), all the p'-local flowing interactions for p' < p are also 
known by injecting the solution for \J\. For example, for p = 2, the bilocal part is obtained 
from ( p9| , p0| ) injecting the solution to (j3"ip . 

To a given order in U one can also perform a loop expansion by expanding the exponen- 
tials of propagators which appear in Q3"T| ) (see Appendix |C.4| for the 0(U 3 ) equation). To 
order U 2 and one loop it reads 

WW) = J?d ii U l {4>)+ fdl' J^Pvi^diPvicf) (34) 

J 

J? = -\jdGl (35) 

Jg, = -\\dG\dG\, (36) 

In order to compute the correlation functions for small wave-vectors q, the strategy of 
the RG consists in performing a perturbative calculation in the theory renormalized up to 
I' = ln(Ao/g). In the favorable cases, the interaction V; flows, from a small initial initial 
interaction V , to "fixed points" in functional space (up to appropriate rescalings) controlled 
by a small parameter (such as the offset from the critical dimension). Once the asymptotic 
large I behaviour of Ui(<f>) is known, one uses the invariance property of W(J) (see (|ll|) and 
Appendix |C5| where this is done in details) to compute the observables. To lowest order in 
U, it is sufficient to keep only the local part in ( |13"D which yields 

W(J) = -J:G:J-[ [ U K e -\K.G^.K e iK.{G:jr + (3?) 
2 Jx J K 

Thus one has for the two-point function (0f</>? )% = 5(q + q')Cfj with 

= G% + J^(K.G q ) i (K.G q )jU l K (38) 
U* = uK e -\K.G r ».K (39) 

and more generally the n-point function (n^2): 

m. ■ ■ ■ & = (n^*^) um m 

with 5 q = (2n) D 5(q). To compute e.g. the two point correlation function at wave-vector 
q, one carries perturbation theory in Uy at a large scale V and it is convenient to choose 
/' = InA/g. To first order in Ui>(4>), one has 

CJ = Gf, + [(KG^lKG^^^e"^^ 1 (41) 

Of course, since W(J) is by construction independent of I the result should not depend on 
the choice of I'. Using the RG flow equation, it can be checked order by order in perturbation 
in powers of Ui that this is the case. 

In the above computation of the two point function, the natural vertex which appear is 
not Ui(4>) but U(4>) = f K e •.^Ui; it is thus interesting to study directly its flow equation. 
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3 Removing of tadpoles and application to 4 



3.1 Modified Polchinski equation 

It is useful for some applications, and in particular to simplify higher orders calculations, 
to get rid of the linear term in the Polchinski equation. This can be achieved exactly by 
introducing the following functional: 

V,(0) = e l^ :G ^V/(0) (42) 
Inserted in (|9|) one finds that it satisfies 

dMffl = \^ Gv ^j^ ■ 9G i ■ ^Wi)U=^(0 2 )U= (43) 

The graphical representation of this equation is drawn in Fig. |2]. Since this equation does not 
contain a linear term, its solution does not contain tadpole-like diagrams. This functional 
has thus several advantages: first it enters directly the computation of any observable, second 
its flow is simpler than the one of the bare vertices. Finally in the context of quantum field 
theory it has the direct meaning of being the normally ordered vertices. 




Figure 2: Graphical representation of the modified Polchinski equation. The point represent 
any vertex, the broken line is a propagator on shell dGi and the full line is a Gi. 



We now perform a multilocal expansion, similar to the one introduced in the previous 
Section, but on the functional Vi(<f>) as: 

V l (<f>) = / Uim+ I Vi(<P x ,<P y ,x-y) + .. (44) 

J x Jxy 

The modified Polchinski equation fl4~3"D can be solved order by order in Ui(<f)). The general 
analysis is performed in the Appendix C.l . Here we give only the result to order U 2 which 
reads: 

dA{4>) = ^/ s e al ^ a2 9 1 9G^ 2 ^(0i)U 1 =^K02)U 2 =^ (45) 

An interesting property of this equation is that it is now local in I. 

Expansion in the number of loops k, restricted to order Uj 2 , is thus straightforward: 

1 +oc 1 

it = Jm k dat = j^dt Jm k+i (47) 
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Note that the multifocal expansion for V and V are not identical, i.e they do not produce 
order by order equations which can be transformed back into each others. However, they 
should yield the same result at the end when calculating universal fixed point quantities. 
This property will be tested and used in the rest of the paper. 

We also give the expression of the bilocal term, as one must check that it effectively 
reaches a fixed point for consistency. It reads as a function of the local term: 

Vi(<k,M = I I {j qx - 1) fdl'e^P&dtfr&MfaA'M (48) 
2 Jx Jo 

3.2 Modified RG equation to one loop and application to 0{n) 
model 

Expanding the exponential of the propagator in (fiS]) yields a loop expansion of the beta 
function to order U 2 . To one loop this gives: 

bUM) = ifdifiWdiMt) ( 49 ) 
I? = \jdGiGl (50) 

As a simple application let us consider the 0(n) model in D = 4 — e with the polynomial 
interaction: 



Ul(4>) = \g2,l^ + ^4,/(0 2 ) 2 + ^6,/(0 2 ) 3 + (51) 

the dimensionless variables being g n j = g n ,i\ 2 • The RG can be performed either using 
Ui(4>) as in the previous section, or in terms of Ui(4>), by which we start. One has: 

f>,(0) = e^r°v 2 W ) = l -a 2 ^ 2 + ^4,,(0 2 ) 2 + -o^) 8 + .. (52) 



and, similarly p(J a n j = a n jA l 2 thus: 



a 2 ,i = d2 t iAf 2 (53) 

a 4 ,z = fl4,/Ap (54) 

a 6:l = a 6j/ A 2 " 2e (55) 

From dijUi{4>) = a 2 ,z5»j + ^-((%0 2 + 20^) + ^f i (^-(0 2 ) 2 + 4</> i i </> 2 ) + .., one obtains the 
RG equations: 

d~a 2)l = 2a v + I D ^{n + 2)~a^~a 2 j + 0(a 2 J (56) 
2 

da^i = ea 4 ,z + I D -(n + 8)5^, + O(a^) (57) 

da 6 ,i = (-2 + 2e)a 6 ,i + 0{a\j) (58) 
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in terms of the single integral I D which has a universal value in D = 4: 

1° = Atf> = At 1 - I C/( A f } C( f } (59) 



2 7, A 2 g2 

= / (2s)-^c'(s)c(s) = -^ + 0(e) (60) 
2 j s >o 4 

with S± = 1/ (8tt 2 ) and here and in the following So is the unit sphere area divided by (2tt) d . 
Thus a^i flows to the fixed value: 

K = , 3 N 16?r 2 e (61) 
4 (n + 8) V ; 

which is universal to this order. We have indicated terms in the above RG equations which 
arise at two loops and yield fixed point for a 2< i and a 6 j with: 

a* 2 = 0(e 2 ) , ~al = 0(e 3 ) (62) 

and more generally a* 2n = 0(e n ) for n > 3. The derivation of this is simple but goes beyond 
this paper Hl6| . This fixed point is unstable in the direction a-i and stable in all other 
directions. The stability eigenvalues to 0(e) read: 

n + 2 . . 

A 2 = 2 e 63 

n + 8 V ; 

A 4 = -e (64) 
A 2n = -2(n - 2) + e(n - 1) n>3 (65) 

The critical manifold of the O(n) model corresponds to a 2 = a* 2 . This corresponds indeed 
to the massless case, since the self energy 

£ 9 = U"(Q) + C(?7 2 ) = A 2 a 2j/ (66) 

vanishes to lowest order on the critical manifold. The instability eigenvalue at the the fixed 
point gives the critical exponent 7: 

7 = |" = l + ^^ e + 0(e 2 ) (67) 
A 2 2(n + 8) 

thus recovering the standard result. One also gets u = — e and r\ = 0(e 2 ). 

Note that to first order in e there is no q dependence of E 9 and no wave function renor- 
malization. This can be incorporated in the method to two loops (see Appendix [A.2| and 

mi 



4 Application to disordered elastic systems 
4.1 FRG equations and universality to one loop 
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4.1.1 The model 



Let us consider an elastic system of internal dimension D embedded in a disordered medium. 
It is described by a single component displacement field u x (x is the internal variable), which 
is either the height function for an interface problem, or the continuous deformation field for 
periodic systems. The energy reads 




.VV2 



-W(x,u x ) + — \u x \ 2 j (68) 
where a short-distance cutoff is implicit. The elastic constant is set to c = 1 here, and the 



mass torn = 0, its effect will be studied in Section jO] . The disordered potential — W(x, u) 
is a random variable which has the following properties (i) W(x,u) = (ii) the potential at 
different x are uncorrelated (iii) the distribution of W(x,u) is translationally invariant in u 
space. Its cumulants read 

W(xi,u 1 ) . . . W(x N , u N f = 5{xi - x 2 ) ■ • • S(x 1 - x^S^iui, . . . , u N ) (69) 

where the symmetric functions S^ N ' satisfy S^ N \ui, . . . , un) = S^ N \ui + u, . . . ,un + u) for 
any u. In particular, the second cumulant is denoted by R(u — u') = S^ 2 \u,u'). In the case 
of an interface, these correlators can be either long-range, e.g. random field, or short-range, 
e.g. random bond. In the periodic case, the cumulants have the periodicity of the lattice. We 
assume parity symmetry S^ N '(—Ui, . . . , —Mat) = S^ N >(ui, . . . ,un). This problem is usually 
studied by introducing n replicas 4> x , a = 1 . . . n, of u x and by averaging over the disorder. 
It yields the action 



5(0) 



Z1 a N>2 ai...aft 



(70) 



Thus the bare S has the general form @ with 



Gl b = ?-c(^)5 ab (71) 
V(0) = / U{4T) (72) 
U{<j>) = ~E R fa ~ ^) - 3^3 E S(3) <h> &) + •■■ (73) 

ab ' abc 

We will consider any cutoff function c(s) such that c(0) = 1 with no loss of generality and 
with c'(0) = for convenience (see below). 

Direct perturbation theory can be performed on this model. One can show that it has a 
well defined T = limit (see Appendix 0). Furthermore, at T = this perturbation theory 
is in fact trivial, i.e the disorder average of any observable is identical to its value in the 
linear random force model, as shown in the Appendix 0. Within the exact RG one can in 
fact escape this well known dimensional reduction phenomenon, since, as we will see below, 
the flowing disorder becomes non analytic. As shown below it yields non trivial results for 
correlations. 
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4.1.2 RG analysis 

We now use the exact RG method introduced above. For now we use the RG equations 
based on the multilocal expansion of V while the other method in terms of V (explained in 
Section |37T] and Appendix |G.1[ ) will be used in Section OL The method with V turns out 



to be more convenient to analyze finite T effects. The / dependence of the Gaussian part is 
implemented by the choice 

(G% = T(B,&= - 2 c(£- 2 )6 ab (74) 
q zi\ t 

where A; = Ae~ l . This choice is particularly convenient here since there is no correction 
to any order to the connected quadratic part (statistical tilt symmetry ||21||). The flowing 
interaction functional V/(0) remains translationally and parity invariant in x space. Since 
translation invariance in the u space is conserved, its local part Ui((f>) remains of the form 
(|75|). In order to obtain fixed points it is convenient to define a rescaled dimensionless 
temperature TJ = TAf~ 2 and rescaled functions 

Ufa) = Ar D W) (75) 

= -iE^(0a-0 6 )-^E^ 3) (^»^^)+--- (76) 
ZI l ab < i .l l abc 

Rtiu) = A?- 4 Ri(u) (77) 
Sl N \u h ...,u N ) = A? N - 2N - D Sl N \ Ul ,...,u N ) (78) 

The general RG equation Q32|) for Ui(4>) implies a set of flow equations for the rescaled 
cumulants (ui, . . . , m^) (since the former is in fact a set of equations for functions 

of a n-dimensional vector <p for any n). The rescalings above have been chosen such that 
these rescaled functions flow to fixed functions denoted R*(u), S^ N '*(ui, . . . , un) independent 
of T. 

An important property of the theory (|70|) is that it admits a well-defined T = limit, at 
least at the perturbative level. This can be seen either by examination of the diagrammatics 
(all negative powers of T in the perturbative calculation of observables are in factor of a 
positive power of n, see Appendix 0), or in the T = dynamics pTSj] . Similarly there exists a 



well-defined T = limit of the set of flow equations for the cumulants. For small e = 4 — D, 
this complicated set of coupled equations can be organized in powers of e. Specifically one 
finds that R* = 0(e) while S^ N > = 0(e N ) for N > 3. This can be seen on the schematic 
structure (ordered in U and T) of the flow equations obeyed by the rescaled Ri, 5$ , which 
can be read off from Appendix |G.4 : 

OU = TU + T 2 U 2 e T + T 3 U 3 e T + ... (79) 



The two first terms reproduce (|31~D since dG oc T, while the third term mimics the 0(U 3 ) 
in the f3 function. Its three U vertices must be linked by at least three propagators because 
of the constraint of locality. Substituting symbolically U = R/T 2 + S/T 3 , where we restrict 
to the two lowest cumulants, one finds 

dR = eR + S\ 2 + R 2 \ 2 + RSe T /T\ 2 + R 3 e T /T\ 2 (80) 
0S = (6- 2D)S + RS\ 3 + R 3 \ 3 + S 2 e T /T\ 3 + R 2 Se T /T\ 3 (81) 
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where we have denoted projections on 2 and 3 replica parts. All terms containing 1/T vanish 
after these projections since a well-defined T = limit exists. We have discarded terms, such 
as diS = TR 2 \3, which (formally) vanish at T = 0. 

One sees immediately on these equations that the fixed R* = O(e) while S* = C(e 3 ). 
This can be generalized by noting that the lowest order (in e) correction to is of the 
form R N \ N thus S^ N > = 0(e N ). To £>(e 2 ) at T = we thus need: 

dR = eR + R 2 \ 2 + S\ 2 + R 3 e T /T\ 2 (82) 
OS = (6 -2D)S + R% (83) 

In this paper we simply perform the 0(e) calculation to which we now turn, for which 
consideration of two replica terms is sufficient. 

We perform the analysis in the T = limit as explained above. The propagator can 

2 

c(3—) ■ A 

be expressed in terms of dimensionless quantities as Gf = T\ J q —Z— e loA i x . At finite T, the 
exponentials of propagators in (|3T|) would reduce to 1 asymptotically at large I. This is also 
true in the T = limit for any I. It is thus a priori unnecessary to include higher number of 
loops within order U 2 . 

Denoting by Mi(<p) the two-replica term contained in the local operator U 

M l {ct>) = - l -Y.R^a-<p b ) (84) 

ab 



the flow equation to one loop (|34|) (using the change of variables V — > / — V yields 

diMi(4>) = (4 - D)M,(0) -\f dl' K v ^2lda b Mi-v(<P)] 2 \2-re P (85) 

1 J ° ab 

where the kernel responsible for the retarded nature of the flow is 

K V = 4J5_„A- ,AP = 2e-^ 1 ' jf d{£ y { l e -*) (86) 

Since c'(u) is typically peaked around u ~ 1 and decreases fast at infinity, one sees on (|86| ) 
that the range of the kernel K/ is also of order one and can be made as small as desired by 
choosing narrow enough cutoff functions. The above RG equation (|85D involves computing 
the contraction: 

^[dadbM^)} 2 = £ [#'(0„ - (j) b f - 2R"(0)R"(<p a - fa)] (87) 

ab ab 

+ J2R"(<Pa-<t>c)R"(<l>c-<t>b) (88) 
abc 

where we have used d a di,M(<p) = 5 a bJ2 c R"{ ( t ) a — 0c) — R"{4>a — 4>b)- The last sum being a 
three replica term, it does not enter the equation for R (it is a correction to S proportional 
to T), which reads: 

miu) = eR^u) + jf' dl' K, Qj#_,,(u) 2 - ^(O)^,,^)) (89) 
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Let us first study the case of periodic elastic systems, with Ri((f>) periodic of period 1. 
Taking the large / limit we find the fixed point equation: 

= eR*(u) + (J^°° dl' K v ) (^R*"(u) 2 - R*"(0)R*"(u)^j (90) 

It is now easy to see that the factor K = J °° dl K; in (p0[), which a priori depends on the 
dimension of space and of the whole arbitrary cutoff function c(s), becomes universal in 
D = 4. Indeed: 

roo rs 

K = 2S D 2~ € / 2 dss~ e c'(s) dtt £ c'(t) (91) 
Jo Jo 

roo 

= 2S 4 dsc'(s){c{s)-l) + 0{e) = S 4 + 0(e) (92) 
J o 

where we used the new variables s = q 2 /2, then t = se~ 21 , kept only the lowest order in e, 
and used c(0) = 1. We denote by So the surface of the unit sphere in D dimensions divided 
by (2vr) D . Thus, to one loop, the FRG equation does not depend on the cutoff procedure. 
It coincides with the fixed point equation obtained |T3], from Wilson's momentum shell 
renormalizat ion . 

In Appendix |CJ.3| , we also mention the result of a two-loop calculation of the beta-function 
in our exact renormalization framework. 

The solution to ([K]) is known to be the 1-periodic function defined by 

R*(u) = -^r (^--u 2 (l-u) 2 ) (93) 



72S 4 V36 

for < u < 1. This fixed point function is non-analytic which is an important and unusual 
feature. It was argued in |14| that this non-analyticity appears at a finite scale. This scale 
Rc — e lc /A can be identified with the Larkin length at which metastability and glassiness 
appears. Taking the fourth derivative at u = of (|89|) yields a closed retarded equation for 
i?f(0) 

dim v (0) = eRT(0) + 3 [ l dl' K^RtAO) 2 (94) 

JO 

In the limit of narrow cutoffs, the equation becomes local and -R™(0) diverges at a finite 
scale. One can show that this feature persits in the non-local equation. 

The case of an interface (i.e a directed polymer for D = 1) corresponds to another fixed 
point where one must rescale the function Ri(u) as follows: 

R\{u) = e 4 %(ue- cl ) (95) 

and we must now determine ( = (9(e) such that ri{v) converges to a fixed point r*(t>). 
Inserting (^) into ([89]) yields: 

dm{v) = (e - 40n(«) + (vrl(v) + ^ dl' K.e" 4 ^ Qr^,(^') 2 - rf_,(0)rf_,(^')) (96) 
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Although the kernel has been modified, this does not affect the results for the fixed point to 
lowest order in e. The fixed point equation reads: 

= ( e - 4C)r» + (vr*'{v) + S 4 Qr*» 2 - r*" (0)r*" '(vj) (97) 

and is thus universal, independent of c(s). This shows that (, which, as shown below is the 
roughness exponent, is universal to one loop. It will be studied below for the random field 
case and in the case of short range disorder it is thus equal to 0(e) to the values given in 



4.1.3 Correlation function 

Let us now compute the two-point correlation function at T = using (|12|) . To lowest order 
in e, it is sufficient to use the first order formula (|4T|). The bare Gaussian part G q vanishes 
at T = 0. We thus get: 

c ab - ~^c(— 2 ) d a wm*eo - A ~ ^ (98) 

where we used that Ri converges to the fixed point R* and small q such that c(g 2 /2A 2 ) = 1. 
In real space, it yields logarithmic growth of the displacements with a universal prefactor 



e 



((u x -u ) 2 ) = —In | Ax | (99) 
18 

In the case of short range disorder (e.g. random bond an for Ising interface) one gets 
instead: 

C ab = -4 = -e A i— ^5 ^T^c ( 10 °) 



This yields to a roughness exponent ((u x — u ) 2 ) ~ |x| 2< > with a nonuniversal amplitude (since 
the FRG fixed point equation ( |9~?D is invariant under r*(v) — > A 4 r*(f/A) and, contrarily to 
the periodic case, nothing here fixed the scale). 

We can now investigate in more details the structure of the asymptotic flow of the various 
higher order interactions (three replica terms and higher, as well as bilocal interaction and 
more). Although this is beyond the scope of this paper, such an analysis is in principle 
necessary for consistency, i.e to ensure the existence of a global fixed point (for all interac- 
tions) and the validity of the result to 0(e). We sketch it here for the periodic case ( = 0, 
generalizations to interfaces being simple. 

We start with estimating the higher cumulants of the renormalized disorder, i.e the higher 
replica components of the local interaction Ui. To lowest order in R, the correction to 
is proportional to Rf and takes the schematic form: 

dS (N) = -{2N + D- DN)S {N) + R!' N (101) 

Graphically, the diagram is made of one loop. We dropped the numerous higher order terms 
in e coming from contractions of various other cumulants than R. One finds that the fixed 
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point S^* takes the following form to lowest order in e 



E s^>(K--M = c NJ} lTt(w ir )- E Rta 1 ---Rta N ) (102) 

ai...a N y ai...aN,b J 

where W a b{4>) = ^abH c R*dc ~ RaL R*al denotes R*"((fi a — (fib), and cn,d is some number 
depending on the cutoff procedure. The last term in the trace has been substracted since 
the product of the N 5's is a N + 1 replica term. For instance, the third cumulant is of order 
C(e 3 ) and reads 

S^*(u u u 2 ,u 3 ) = c 3<D Sjm UuU2tU3 [R*"(u 2 - u 3 )R*"(u 3 - Ul )R*"( Ul - u 2 ) (103) 
-3R*"{ Ul - u 2 )R*"{u 1 - u 3 f + 3R*"{0)R*"{u 2 - ujR*"^ - m)] (104) 



where c 3)4 is computed in [|nj and reads: 



^-IjT*^ <«■> 



Now we check that the bilocal part has a well defined fixed point. Its expression is given 
by (|29"D, where at T = 0, the exponentials should be expanded at most to first order. The 
zero-th order term yields three replica terms, while the first order term yields two replica 
terms (as well as a correction proportional to T to three replica terms which we can discard 
at T = 0). Thus we get 



vmum = -§- 2 Y.vr\j-) - f^z vr ] (.j- > doe) 

Z ' 1 ab 1Vl UJ abc IVl 



where we have explicitly separated two and three replica terms respectively: 
V^ 2 \q) = ~ J k ( k2{k \ ~ )2 (c(k 2 /2) - l)(c((fc + q) 2 /2) - 1) - idem q=0) 

(E R abKb ~ R "(°) z2( R 'ib + Kb)) 
ab ab 

--( f ±- 2 c(k 2 /2)) f da(- - l)c'(aq 2 /2)(2 £ &*&L + 



il r>///2 \ . 

2rep 

abc,a^baj^c 



V abc 

where stands for i?'(0 l a — 01^) etc.. The bilocal term thus has a scale invariant fixed 
form of order e 2 and is a well-defined function of q = q/Ai with no divergences. 

More generally, we conjecture that there is a fixed asymptotic form for all multilocal 
interactions V^' which can be explicitly written as a sum of properly rescaled multi-replica 
terms as 

V^\^.4 P ^..x p ) = A? p eC^ E V^ c \{ ( j )a ^ : ;;lki^.A lXp ) (107) 

c>2 C - 1 ai,..o c 
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where the number of replicas c corresponds graphically to the number of connected compo- 
nents. The consistency of the method demands that the V\ flow to well defined fixed points, 
perturbative in e. It is indeed natural to conjecture that in this theory there is no wave 
function renormalization. 

We can now come back to the calculation of the correlation function. Although for 
convenience we have computed C q ab from the theory with I = ln(g/A), this is unnecessary. 
As discussed above, the existence of a fixed point with well-defined functions of q = q/Ai 
implies that the e expansion of C q ab at fixed q <C A should be of the form: 

C ab = ^AfE(ff) (108) 

where the dimensionless self energy S(g) depends only on q and e. Since C 9 ab is I independent, 
it implies, taking the derivative, that S(g) = C e q e and thus 

C ab = C e q~ D (109) 
C e = 2 -^e (110) 



to lowest order in e, as in [23|. The form (|109|) should be valid to all orders in e if the 



hypothesis about the fixed point formulated above are satisfied. 

4.2 FRG at finite temperature 
4.2.1 Renormalization equations 

Since in (|3T| ) the terms in the exponentials containing the temperature go to zero as Af ~ 2 
one can first study the effect of temperature, compared to T = by looking at the linear 
term. Up to order TR and R 2 (i.e to one loop) the RG equation thus reads: 



dMu) = eR^n) + m'liu) + f dl' K v {\r!U,{uY - R'i v {^Lv{< 



with 



T, 



T I Afd< q 



2TS A A 2 j c{s) + 0{e 

Js>0 



2A 2 ' 



(111) 

(112) 
(113) 



In the case of sharp cutoff this equation has been studied in p5j [26|1 . It was found that 
at fixed u, Ri(u) converges to R*(u) but that temperature rounds the cusp of the T = 
solution in a boundary layer of size u ~ TJ. As in [26| we look for a solution of the form: 



ir 



T? 



e 2 x 2- - ^ 



(114) 
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Here, x (°f order e°) is defined by ex = R*"'(0 + ) and thus H*""(0) = —1. One has also 
Hi(0) = H"(0) = 0. Injecting (|114|) into ( [Lllj) and identifying the leading order in 2], one 
gets: 

— = if*"(x) + - I K ie il H*'\xe- 21 ) 2 (115) 
2 2 7/>o 

with J i>0 Ki = 1. This equation can be solved iteratively in n for H*(x) = Y, n >o a n% 2n / (2n)!. 
One has a x = 1, a 2 — — 3 but higher a n 's are non-universal. The large x behavior of H* is 
universal and given by H*(x) ~ x. In the case of sharp cutoff one recovers fl2"E 



iT"(x) = Vl + x 1 - 1 (116) 

This result should be further examined by consideration of consistency within higher loop 
corrections, which goes beyond this paper. 

The most important result of the section is that the following relation between the finite 
temperature solution and the T = solution: 

lim TiiC'(O) = C7?;'"(0 + ) 2 (117) 

(— >+oo 

where C = 5 4 holds irrespective of the cutoff function, and thus is determined by the T = 
fixed point. This property will be used below. 

4.2.2 Calculation of universal susceptiblity fluctuations 

It was noted recently that a signature of glassy behaviour in a disordered system was the 



large sample to sample fluctuations of the response to external perturbations p7| , [29| , |3 
These are described by the following suceptibility: 



X 



~f ({d a u x d aU y) - {d a u*){d aU y)) (us) 

1 Lj J xy 



in a finite system of size L, which measures the response in a given sample to a field coupling 
to Vm (e.g. the tilt or compression a flux lattice, or the compression response). The (X) 
denotes the thermal averages in a given sample. These have been studied in connection with 
mesoscopic behaviour of disordered systems |HJ. Here we have considered only the trace 
of the response tensor (extension being straightforward). To perform the calculation in the 
replicated theory we define: 

C$ = ~ fjd a «° a dffl® (119) 

C abcd = ^-4 / (V< ■ Vu\Vu z c ■ V?4> (120) 

1 Li J xyzt 

We now compute respectively the first and second moment of the sample to sample fluctua- 
tions of the susceptibility. They read: 

X = C aa -C ab = l (121) 

X 2 — C a abb + Cabcd ~ 2C aa bc (122) 
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where X denote disorder averages and a, b, c, d take values all distinct from each other. Note 
the well known property |21| that the average susceptibility is identical to the susceptibility of 
the pure system. We now compute C a b c d to lowest order in e. Only zero and one loop graphs 
involving respectively R" (0) and R""(0) contribute. Interestingly, due to the quadratic nature 
of the term proportional to -R"(0) the zero loop graphs cancel in x 2 -, as can be easily seen 
since C a b C( i = C 2 b + 2 Y^osi^^f) 2 f° r an y gaussian theory (performing the Wick contractions). 
One is left with: 



Cabcd = (n5 abcd - (5 bcd + 5 c da + 5 d ab + 5 a bc) + (5ab5 c d + 5 ac 5bd + 5 ad 5 bc ))AR'i'{Qi){l22) 
where A = L~ 2D J w (J xy VG X ~ W ■ VG y ' w ) 2 and thus 

(Axp =^-x 2 = AR';"(0) ~ R'r(0)CL 4 - D (124) 

for a system of finite size L (see also |SIJ for a similar result in straight perturbation theory). 
Note that one can equivalently study the perturbation of an infinite system (i.e L — > oo first) 
by a periodic external field of wavevector q ext . In that case A = q^ A . Thanks to the exact 
RG equations at finite T and substituting I = lnl we obtain the mesoscopic susceptibility 
fluctuations at low temperature as: 



(A X ) 2 = C— (125) 

where 9 = D — 2 + 2£ is the energy fuctuation exponent and C = 0(e 2 ) for a periodic system 
(£ = 0) and C = 0(e 4 / 3 )a 2 / 3 for an interface in random field disorder (see Section |4.3|) . 
This result, derived here through exact FRG calculation, is consistent with the droplet 



picture |32]| . Indeed the second moment of the susceptibility fluctuations is dominated by 
the rare configurations of disorder (of probability pdeg <^ 1) with two almost degenerate (i.e 
within 0(T) in energy) ground states as follows: 



(A X ) 2 ~ p dcg (5 X ty P y (126) 



where pd eg ~ T/L 9 and the typical fluctuation is Sxt JP ~ T l L D L 2D 2+2< * from |118| . One 
thus recovers the above result since 9 = D — 2 + 2(. 



4.3 Interface in a biased random field and toy model 

In this Section we study the model ( |68D in the presence of a mass term m > 0, which confines 
the fluctuations of the displacement u x . We consider two cases (i) random field disorder (ii) 
periodic disorder. A physical realization of (i) consists in a domain wall separating the ± 
phases in a ferromagnet, submitted to a random magnetic field. The magnetic energy of the 
interface, assumed without overhangs, is: 

E(u) = 2 J d D x£ du'h(u',x) (127) 

Thus the effect of the mass term corresponds to applying an additional field gradient 
h(u,x) — > h(u,x) + m 2 u/2. Note that this field gradient can either stabilize (m 2 > 0) 
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or destabilize (m 2 < 0) the domain wall. We will study the approach to the critical value 
m 2 — > + . The case (ii) is of interest when studying the competition between disorder and 
e.g. a periodic potential. In the phase where the periodic potential is relevant it is natural 
to approximate it by replacing it by an harmonic well (see e.g. [|33[|). 

Here we examine only ground state properties (zero temperature). We show that the 
disorder induced fluctuations of the displacement u x is described, as m — > by a universal 
scaling function of the form: 



u i u -i = m~ a F[cq 2 /m 2 ] (128) 

which we determine to lowest order in e = 4 — D. Note that c can be measured from the 
thermal connected correlation, which is unchanged by disorder for models like (|68"D which 
possess the statistical tilt symmetry. 

4.3.1 RG equations in presence of a mass 

In this Section it is more convenient to use the RG equation resulting from the multilocal 
expansion on V, which is local in I to this order. Since we are studying T = we set Ri = Ri 
in the following. The RG equation reads: 

diRi(u) = Ji {^R'l{uf - J2f(0)Jtf(u)) (129) 
Ji = -2 / dG q {Gl (130) 

Jq 

Using the action Si, the T = correlation function reads to lowest order in R: 

;*Fi =_«.(„) (£<^>y (131) 

\cq 2 + m 2 J 

It is easy to transform ( |129j ) into the RG equation in the absence of a mass. Using the 
change of variable: 

R l (u) = c 2 Al^R m (uA c m ) (132) 
one finds that Ri satisfies the m = flow equation: 



d t R t {x) = (e - K)Rt{x) + (xR' t {x) + S D ( \r"{x) 2 - R'l{Q)R'l{x) ) (133) 



2 



provided that t{l) satisfies 



2 



which is integrated into: 

e^-1 r°° D _ Xj c(q 2 /2) 2 - c(q 2 e 2l /2) 2 
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with a = The function t(l) is increasing and bounded. Its limit t(+oo) 
is given by: 

e etoo _ l f oo n _ c{q 2 /2) 2 



ii {q z + ay 



for D < 4 



(136) 



and it diverges for m — > as: 



e et °° ~ e / g D -^g 7T - — ^ 
Jo (q 2 + a) 2 



a-tfl-|) 2 



2 y sinf 



In D = 4 it diverges as: 



We now distinguish the two cases. 



4.3.2 Random field 



1, ,1, 

too ~ - ln - 

2 a 



(137) 



(138) 



In that case the correlations of the potential are (W(r, u) — W(r',u')) 2 = —25 D (r — r')R{u — 
u') with (from ( |127| ), R(u) ~ — a\u\ at large u. In the massless case it is known that the 
FRG to one loop reproduces the purely dimensional result uqu t ~ <j 2 /3 c -4/3 r 2(4-D)/3 with a 
roughness exponent ( = (4 — D)/3. From this we expect, in the massive case, the small m 
behaviour: 



It is known that one must fix £ = e/3 to obtain a reasonable fixed point. From the 
above equation, the (reduced) correlator of the force A t (x) = —R'l{x)Sr> then satisfies 

dtMx) = 6 - (xA t (x))' - l(At(x) - A t (0)f 
and flows for t — > oo, to a fixed point A 0O (a;) given in terms of a function [III] implicitly 
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Aoo(0) 

2 



A oo (0)i/(x V 3 Aoo(0); 

fe) 1/3 (/- + ^Ao(^ 2/3 
y(x) — 1 — lny(x) 



with 7 = Jo 1 Vy- 1 -lny = 0.5482228893.. Note that f+™ dx A t (x) is t independent and 
thus equal to /+£ cftc A (x) = -f£ / R' ' = 2a^-. 

Putting this together with (|131|), (|137|) yields the result: 



(139) 



= a 2 l z m- D ~^c D l & F D [cq 2 /m 2 



F D [x] = C D 
C D 



(1 + x) 2 
/ (4tt) d / 2 \ 1/3 

Ur(f)7 2 J 



+ h.o.t 
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with ( = e/3. Note that the universal scaling function must behave at large x as Fd[x] ~ 
x^~ D / 2 } and is determined here only to order in e. 
From this one also finds the local fluctuation: 

which is also universal. The fact that this quantity is dominated by large scale fluctuations 
can be seen from the convergence of the integral (£ > 0) . 

The calculation can also be performed exactly in D = 4. One finds: 

/ 1 x -1/3 /A~2\ 1/3 

A^ = A 2 Vln- 1^- (141) 

V mj \3^J (1 + ^) 2 1 ' 

The values we find for Co, which vanishes as Cd ~ 3.5246 e 2//3 as D — > 4~ are as follows 
C 3 « 2.40653, C 2 « 1.91006, Ci w 1.30416, C « 0.82157, remarkably close to the exact 



result in D = [H C = 1.05423856519. 



It is also useful to compare these results with the Gaussian variational method with 



replica symmetry breaking. Extending the calculation of |27| to the non zero mass case 



which is done in Appendix [E], we find the same form as (|139|) with 



F M -r'\ 1 I 6{2 ' e) 1 / + °° dy V ~\ (142) 
1 9 

C' D /C D = (— 7 2 ) 1/3 = 1.04708 (143) 

7T 

and 9 = (6 — e)/3 and ( — e/3. Since as e — > for fixed x the second integral is subdominant, 
the leading order in e are identical and the amplitude of the RSB solution compared to the 
FRG solution is C' d /Cd as e — > 0. 

4.3.3 Periodic case 

In the case of a periodic system with period a, one gets a fixed point function with ( = 
which reads: 



e ( a? 



It yields: 



and in Z) = 4: 



A 00 (u) = -l--u(a-w)J (144) 



„ n/9 n (47r) D / 2 1 , . 

u iu~i = a 2 c D/2 mr D K , — - 145 

36r [f] (i + C S) D/2 

„ ^-a 2 H-) (146) 
36 m 



27T 2 2 m 



u i u -q = aV- — , — (147) 

9 1^(1 + ^) 2 
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5 towards two loop FRG 



The exact RG method allows to compute quantities beyond the lowest order in e. It can be 
carried either at T — > for fixed system size (T = limit) or at finite T. Solving the exact 
RG equation at T = requires to follow non analytic functions. This is a difficult question, 
e.g. distinguishing the various cumulants in the local part demands a special procedure that 
we have developed. This is discussed in |]TBl ITS . 



These problems do not arise at T > where the singularity is smoothed within a bound- 
ary layer (at one loop see Section fj~2[) . In the Appendix [F| we have used the exact RG flow 
to third order (given in |G.4| ) and obtain the two loop exact FRG equation for the second 
cumulant Ri(u) at T > 0. At large I the effective temperature TJ — > and one recovers 
an "effective" zero temperature equation. This equation differs from the one obtained in 
44| as it contains a new "anomalous" term of the form XR"'(0 + ) 2 R"(u). We find that the 
coefficient of this term is universal with A = 1/2. Interestingly this value is consistent with 

The question of calculation of correlations and 
Let us mention here that our exact 



the renormalizability arguments given in [[15 



of their universality to two loops is rather delicate |TB 



-c 



RG fixed point equation depends only on one non universal coefficient K , which vanishes 
in the case of sharp cutoff. The solutions of this equation, given in Appendix |F|, exhibit the 
property of a non zero value of R'(0 + ), referred to as a supercusp since it is a stronger non 
analyticity than the one loop one (R"'(0 + ) ^ 0). This feature is unpleasant as it naively 
yields (by perturbative expansion) additional divergences as T -1 / 2 , a sign of possible frac- 
tional dependence in e (a related phenomenon is discussed in [0]). On the other hand since 
its magnitude is proportional to K it is absent in sharp cutoff calculations. This problem 
and its relation to the structure of the boundary layers at high orders is further examined 
in IITB . 



6 Conclusion 

In this paper we have introduced a systematic method which turns the exact, though ab- 
stract, RG functional equation of Wilson-Pol chinski into a tool for concrete perturbative 
calculations to any number of loops using arbitrary cutoff functions. The strategy was to 
explicitly integrate out all non local interactions, which can be expressed in terms of the 
local part alone, order by order in the local part. In the process we have preserved the 
exactness and the controlled nature of the original Wilson-Polchinski equation. Indeed, no 
approximation was made, and the resulting RG equation for the local part, as well as the 
expressions for the nonlocal ones and for the correlation functions, are formally exact order 
by order in an expansion in the local part. This expansion will be useful for theories where 
the local part is small, i.e when it is controlled by a small parameter (e.g. the shift from 
the upper critical dimension) and when the RG equation admits a perturbative fixed point 
in this parameter. We have considered here theories with a bare local interaction and a 
fixed point for the local part, e.g. as in the 0(n) model, but the method is more general 
and can be extended to theories where the bilocal part of the interaction serves as the small 



parameter, e.g. for self avoiding manifolds j34|. In a sense, the exact RG in the operational 



form presented here directly translates the ideas of Wilson and provides explicit checks of 
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universality. 

In addition to presenting the method formally to all orders, we have derived the explicit 
RG equation for the local part up to third order. Further expanding in the number of loops, 
we have explicitly given the coefficients up to two loops and third order. Two distinct, 
although equivalent, methods have been presented, depending on whether one considers 
the Wick ordered functional or not. Each method has its advantages: the Wick ordered 
method yields apparently simpler (less nonlocal) RG equations, but it is not always the 
most adequate (e.g. for the finite T one loop FRG analysis). Although the present paper 
contains all the material necessary for two loop applications (e.g. for 0(n) and for the 
FRG) we have preferred to defer giving the detailed calculations and results to a companion 
paper |Tj|. In particular we have sketched here the simple extensions needed to deal with 
the so-called wave function renormalization which arises in e.g. the 0(n) model to two loops, 



examined in more details in HI6 



We have thus considered here mostly one loop applications. The first one was a simple 
check to recover the one loop exponents of the 0(n) model. The second application was to 
the theory which describes elastic systems in random potentials. It was previously analyzed 
through simpler Wilson momentum shell integration H] but the rather unusual nature of 
the theory (n — > limit, non-analyticity) made it important to verify explicitly that the 
results are universal. Also universality in disordered systems is rather less established than 
in pure systems, especially in the T = limit where it is known to fail in some cases. Thus 
we first derived the T = one loop RG equation for the second cumulant function R(u) in 
arbitrary cutoff scheme and found that its coefficients are universal to this order. This yields 
the universality to 0(e) of the roughness exponent ( of pinned interfaces. In the periodic 
case, we also explicitly verified that the correlation function contains a universal amplitude. 
Similarly, we computed the scaling function of the ground state deformations of a confined 
interface in a random field and found a universal result. This quantity can be experimentally 
measured in disordered magnets in the presence of a small additional field gradient. 

Although temperature is formally irrelevant (the dimensionless temperature flows to zero) 
it is well known to be "dangerously" so. Our exact FRG at T > shows that although the 
"boundary layer" , i.e the detailed asymptotic form of the cumulant R(u) for u ~ T} is nonuni- 
versal, some of its features are universal, and in particular we were able to extract from it the 
universal divergence of the mesoscopic fluctuations of the suceptibility A%. The divergence 
of this quantity, which is dominated by rare almost degenerate low energy configurations, is 
an accepted unambiguous measure of "glassiness" in a disordered system and is measured 
in experiments, e.g. in microsize vortex systems. 

Some of the peculiar features of the theory of pinned elastic systems have been also 
discussed. We have found it useful to give a detailed diagrammatic proof of the triviality of 
naive perturbation theory, as we have not seen it explicitly in the litterature (though more 
general statements about dimensional reduction appear in a number of other works). We 
have discussed how the non analytic nature of the theory yields non trivial results. 

Pinning of disordered media thus provided us here with one example of a problem where 
exact renormalization is needed to get insight, as no field theoretical description is yet 
available. The reason for it is that one must follow in principle a rather complicated object, 
the full probability distributions or the disorder, or equivalently the whole series of cumulants. 
The method seems thus promising for other problems with similar features, such as random 
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Sine Gordon models [p5| ]. It is interesting to note that while presenting the FRG method, 
Fisher pointed out (Ref. 12 in |14j]) that the momemtum-shell RG "suffers from pathologies 
due to the sharp cutoff", and that the cusp in R{u) "requires a careful analysis of the full 
renormalization group" . Through the use of the exact RG method presented here, we provide 
a simple way to integrate explicitely and exactly what is left aside in the traditional RG. 
We are thus able to control the approximations of former approaches and to perform new 
calculations. Furthermore, thanks to our general framework expressed in terms of any cutoff 
function, the universality of the results is checked. 

Let us close by noting that the application of the multilocal solution to the exact RG 
equation seems promising also to study other disordered problems, or even give a new per- 
spective on simpler pure problems. For instance one could apply it to wetting problems 
taking into account the nonlinear part, or to the roughening problems to improve on pre- 
vious analysis using uncontrolled projections methods j36| The multilocal expansion allows 
also interesting extensions to theories with bilocal bare action, such as polymers, mutually 
interacting or with disorder. Finally, it is also worth studying more closely the set of exact 
solutions to the Polchinski equation presented in this paper (Appendix ||). Some of these 
extensions will be explored in future publication. 



A Invariance properties of generating functional and 
renormalization equation 

In this Appendix we give a concise derivation of the exact invariance properties of the gener- 
ating functional of correlation functions under coarse graining. These properties provide the 
basis for developing exact renormalisation procedures of the Polchinski type. In the second 
part we generalize the framework to include additional field transformations, such as rescal- 
ing. This extended framework is suitable for theories where wave function renormalization 
must be included (see |]TE|]). 



A.l Invariance under coarse graining 



We use only the two following properties of Gaussian averages. The notations are the same 
as in the body of the paper. First, transformation under a change of variable <p —>■ <p + ijj for 
any field ip in the functional integration over <p yields: 



[Am 



G 



e 2 



■tfi-.G- 1 :^ 



+ V0 



for any functional A{4>). We will also use the composition property: 



Wi + ^)] G J G = W)U 



-G 2 



:i48i 



(149) 



where in the l.h.s. the average over 0j is performed using Gaussian correlations Gj. 

Using successively the shifts 0i — > 0i — G 2 : J and <p 2 — ► (P2+G2 '■ J yields the fundamental 
relation 



Go 



3 J:(01+4. 2 )-V(4>i+</> 2 ) 



G 2 



(150) 
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-\J:(G2+G2.G- 1 :G 2 )J 



3 J:(l+G 2 :G ] - 1 ): 9 ii 



-V{<t>i+<h) 



G 2 



(151) 



where we denoted G2 = Gq — G\. Thus, if V\ is coarse-grained transformed of the interaction 
V, defined by 



-V(</.i+0 2 ) 



then one has 



J-d>-V{4>) 



,2 J:(Gq— Gq:G 1 :Gq):J 



G 



Gq—Gi 



,J:G :G^:4>l p -Vi(<t>l+4>2) 



Gi 



(152) 



(153) 



We now use this property of Gaussian integrals as follows. One defines a family of actions 
Sg{4>) and their associated generating functional Wg(J) 



S G (<P) = -(f) : G^ 1 : + V G (0) W g (J) = In 



J:<f>-V G (<f>) 



G 



(154) 



They are indexed by the matrix G and we choose them to be related by the coarse graining 
operation ( |152| ) where G plays the role of Gi, namely: 



-V Go (0+V2) 



Go-G 



or, equivalently in a differential form, the Vg satisfy the "RG equation" : 

5G 



(155) 



(156) 



obtained by differentiating (|155|) with respect to G. The coarse-graining equation (|156|) , 
read along a given path 1 1— > Gi, is the Polchinski equation in its "diffusive" form (|r4"|). 
It is easy to see from (|153|) that this choice of a family Vg implies the property: 



Wg(J ■ G : G~ ) independent of G 



(157) 



where we have defined the interaction part Wg(J) = Wg(J) — \ J '■ G : J . It allows to relate 
correlations within any member of the family Sq, i.e under coarse graining. 



A. 2 Generalization including rescaling and change in Gaussian 
part 

The previous properties can be extended to a larger set of transformations which include 
simultaneous (i) coarse graining (ii) linear transformation of the fields (iii) redefinition of 
the Gaussian part (such as needed to absorb its possible renormalization) . It is based on the 
following properties of Gaussian integrals. One defines: 

W Gy {J) = \n[e J --f- v ^} G (158) 

The first property correspond to performing an arbitrary linear transformation on the field: 

W G y{J) = W M ^ G :M-\V(M:d J ' M ) ( 159 ) 



2$ 



valid for any G, V, J, M. The second property is simply the identity obtained when redis- 
tributing the Gaussian part: 

W G>V (J) = ~\tt ln(l + H : G) + W {G - 1+H) - 1>VW _ J) (160) 

valid for any G, V, J, H. 

First, let V G (0) satisfy the RG equation ( |156[ ). Then from the previous Section we know 
that W G y G (J : G' : G~ l ) - \J : G' : G- 1 : G' : J is independent of G for any J, G'. Setting 
G = Gi and G' = Go one gets the Polchinski equation and one can compute W(J). One 
now defines 

V g ,m(0) = V G (M : <f>) + l -<\> : (M : G" 1 : M - G" 1 ) : <j) (161) 

Using the above properties one has 

W G y G M {J : G' : G- 1 : M) - -J : G' : G" 1 : G' : J - ^Tr ln(M _1 : G" 1 : M -1 : G) (162) 

is independent of G and M, for any G', J. 

We have used the two invariances choosing H = M" 1 : G" 1 : M _1 - G" 1 leading to the 
intermediate formula: 

Wo,v(^) = —TcHM- 1 ■ G- 1 : M- 1 : G) + ^v^j+i^-^-o-i)^ : M) (163) 
Defining a new family of functional indexed by I as: 

V/ = V G; , M; (164) 

and symmetric matrices G, M one finds that the functional Vi now satisfies a new RG 
equation 

diV = -7— : M _1 : «9M : + : G" 1 : M -1 : dM : (165) 

(70 

1 <9 2 V 19V <9V 

_i T r(<9G: + :«9G: ^- (166) 



which contains additional terms. In particular the quadratic piece can be used to absorb 
"wave function renormalization" terms, so as to keep V; small. Once this equation is solved 
the correlations can be related within any of the corresponding Si theories using the above 
invariance property (|162|) of Wg,v g m (J)- This will be further exploited in |16 . 



B General properties and exact solutions of Polchinski 
equation 

Let us first mention a few general properties of (§,0). For the class of cutoff functions 
(H) used in practice, the diffusion tensor in fll4]) is positive c'(s) < (but not definite 
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There are some exactly formally con- 
Since (0) is a diffusion equation, 



positive since there exists modes with d\G\ = 0] 
served quantities, such as f^e~ Vl ^' and e~ Vl ^ ^ 

it satisfies a H-theorem of increase of the "entropy" Si = L Vi((J))e~ Vl ^\ which flows as 
diSi = -§ j^e- v '^ : diGi : ^e _v, W > and is compatible with the fact that RG trajec- 
tories do not have limit cycles. Finally, since QI3]) is a linear equation, if we now a set of 
solutions Vf (0), then any superposition such as 

V,(0) = -lnJ^cVW (167) 

is also solution. 

This can now be used to construct non trivial exact solutions to the Polchinski equation. 

The simplest family of exact solutions is of course the quadratic potential, for which one 
finds the solutions 

Vi{4>) = ~(0 - V) : M l : (0 - V) - ^Tr InM, (168) 

M { = (Mo" 1 + G - d)' 1 (169) 

where ip is an /-independent field, with (<j>) s = (1 + M" 1 : G -1 )^ : V- 

A much less trivial family of exact solutions of Polchinski equation is obtained by super- 
position of gaussians, i.e of quadratic potentials. It reads: 

Vi{<f>) = - In Y.° a e-5(^ Q ) :M r : (*-^)+5 TrlnM r (170) 

a 

with arbitrary constant coefficients c a and each Mf satisfies ( |168| ). This is somewhat rem- 
iniscent of a decomposition into "pure states" and is clearly of interest to describe low 
temperature states in pure models (in phases with broken symmetry) or in disordered mod- 
els and glasses (with many metastable states). It is an interesting question to ask, quite 
generally, whether this family of solutions can in some cases be an attractive manifold in a 
larger functional space, or whether one can carry perturbation around this subspace. These 
and related issues will be discussed in a future publication ]45| 

A generic property of these solutions to Polchinski equation is to generate cusp singular- 
ities separating the "pure states" . This can be seen directly above since negative curvatures 
tend to increase in absolute value (see (|168|) and is presumably a very general mechanism. 
It can also be seen on the simple example of the zero- dimensional toy model. There, the 
field is a real number and 

Z = T d<Pe-^- v ^ (171) 



where V is an arbitrary function. One can introduce Gi = G — I, and Vj(0) for this model, 
verifying 

W) = 5 (tf (0 - W) 2 ) (172) 

with initial condition Vq(0) = V ((/>), and one can integrate up to I = G. One has e~ Vl ^ = 
-\ [<>-<•! rp ne evo j u tj on f Vi((f)) is that the curvature Mi = VJ"(0), which would obey 
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( |168| ) Mi = Mq/(1 + MqI) for a quadratic hill or well, diverges at a finite I for maxima and 



decreases as 1// for minima of V{4>). Thus the landscape Vi(4>) develops cusps, encoding for 
discontinuities in the force — V'(0). 

In the case of a periodic landscape, the natural superposition of gaussian solutions is the 
Villain potential V(0) = — ln^ n ce~^ n ^ ^ 2l \ In these Sine gordon type potential, as well as 
in the 2D XY model, it is a well known property that the renormalized potential converges 



towards the Villain form at low temperature as found in [37] from the Migdal Kadanoff RG 



see also more recently P5|). The detailed behaviour of the RG flow can be studied in a 



more controlled way using the method presented of this paper [45 



C Multilocal expansion and higher order RG equation 

In this Appendix we derive the systematic multilocal expansion and obtain the RG equation 
to higher orders. We give a detailed presentation for the functional Vi(<f>), which is simpler, 
and give explicitly the corresponding RG equation to order Vf and up to two loops. Then 
we simply sketch the result for the same procedure applied to the functional Vz(0), which is 
more involved and will be presented in ||16| . 

C.l Multilocal expansion for V 

The tadpole- free functional V;(0), defined in d42|) , can be written as a sum of multilocal 
interactions 

H(0) = E / V<*\<j> Xl ..<l> Xp ,x 1 ..x p ) (173) 

p>0 Jxi..x p 

Note that we are not even assuming here translational invariance. The translationally invari- 
ant case discussed in ( |2.2|) can be recovered by setting V l (0i, x%) = Ui(4>i), V l (0i, 02, Xi, x 2 ) 
V;(0i,02,£i — x 2 ) etc.. Since we want to impose that each p > 1, has zero local part 
(this is sufficient for our purpose), we define (extending fll9| , [20D ) respectively the projection 
operator Pi which projects a p-local interaction on a local one, and the projection operator 
P± which transform a p-local interaction into another p-local interaction as: 

(PiA)(cf>,t) = f 5(t- Xl + --- + Xp )A(<p,.4,Xi,..x p ) (174) 

Jx\..x v p 

(Pi A) (0i, ..0 p ,xi, ..x p ) = S(xi - x 2 ) ■ ..S(xi - x p ) (175) 



x f 6{xi- Vl + " ' + Vp )A(<f>i . . . p , yi . . . y p ) 
J yi---y P V 



The property 



X\,..X. 



{P 1 A){<p xi ,..<f) Xp ,xi,..Xp) = J{PiA){<f) t ,t) = J t A{<f) t ,..<f) t ,t,..t) (176) 



ensures that one can choose the V^ p \ p > 1 in the decomposition ( |173j ) to have no local part, 
i.e: 

PiV^ = PiV^ (p) = (177) 
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for any I by applying Pi and 1 — P\ act on both sides of the Polchinski equation. 

Since the modified Polchinski equation (|43|) concatenates two operators, it is then easy 
to see that if the satisfy the following set of equations: 

d l vV{<t>,t) =4EE/ S(t - Xl + • • • + x * ) ( ?"« G ' a " 1 ---'Bh«dG l d* 1 "* 

2 p>0q=l Jx l" x P P 

V®(<j> xi . ..<j> Xq ,x 1 . ..x q )V^\<p Xq+1 . ..<p Xp ,x q+1 . ..x p % i=4> (178) 

$v<«« ...,,,,.,■,...,,) is(i - Pi) g e ^-'^' +l -^i... 99Gi r 1 -^ 

2 9=1 

Vb\<f> Sl . . . <j) Xq , x x . . . x q )V^- q \<p Xq+1 ...<p Xp , x q+1 ...x p )iorp>l (179) 

then ([43D is obeyed by Vj(0). Since we prefer to work with symmetric functions we have 
defined the symmetrization operator: 

SP (p) (0i...0p,Xi...X p ) = -J ^2 B(P \4>*(1) ■ ■ -(pa(p),X a (l) . . .X a{1) ) (180) 

we have also defined the following shorthand notations: 

d x...q dGdq +i... P = ]T ]T 9Gl aX0 d?df (181) 

a=l /3=q+l 

It is easy to see that if V^ 1 ' is considered formally as "small" in some sense (e.g. controlled 
by a small parameter such as e) then one can integrate exactly these equations order by order 
in and check that = 0(V^ P ). More precisely, to a given order one can exactly 
integrate the equations for higher point functions and reduce to a single equation for 
This is the procedure that we now follow. The structure to the lowest order 0(V^ ) is 
simply a closed equation for of the schematic form: 

= PT(V (1) * V {1) ) + 0(V {1)3 ) (182) 

To next order OiV^ 3 ) one needs to solve the coupled set: 

d t V^ = P[iy {l) * V w + V w * V {2) ) + 0(y (1)4 ) (183) 
d t V^ = (1 - Pi)(F (1) * y (1) ) (184) 

The second equation is explicitly integrated which yields [V^] which is then substituted 
in the first equation, producing a closed equation for This procedure can be extended 
to any order in We now give the explicit calculation. 

C.2 RG equation up to order V^ 1 ' 

To (lowest) order V^ 1 -* 2 , the beta function is local in I as the modified Polchinski equation 
itself and reads 

d^ (1) (<M) = \ l ^-^^)e alGi °VdG^V%i,xi)^ 
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up to terms of order 0(V^ 3 ). 

To next order V^ 3 , as explained above one first compute the bi-local operator as a 
function of V W . Its flow equation to the necessary order reads: 

fyV^fafaxM) = ^e olGi92 9 1 9G / 9V«(0 1 ,x 1 )\/«(0 2 ,x 2 ) (186) 

-8{x x - x 2 ) I 8{ Xl - yi±yi)e 9lG ' a2 d 1 dG l d 2 V^( ( p 1 ,y 1 )V^(4>2,y2) 
JyiV2 * 



up to 0(V^ 3 ) terms. Integrating <9 M V^ 2 )(0i02, X1X2) using ( |186|) from to I and substituting 



the result into the equation for one finds the RG equation of the local part of the 
interaction to order : 

dtVMfat) = \l 5{t-^^)e dlG ^d l dG l d 2 V^\Kx l )V^\^x 2 ) (187) 



2 



Xl'X-2 



2, J X1X2X3 o \J 



, '' Jo ^s 1 ^ M a B yi 1 H^,a:i)U^ 1 )(0 2) x a )|^- < y( a ; 1 -x a ) / s( Xl -^^: 

J yiy2 ^ 



e 

up to C(V (1)4 ) terms. 

C.3 Translation invariant theory and loop expansion 

In a spatially translational invariant theory the local interaction does not depend explic- 
itly on the space variable t, V^'((f),t) = Ui(<f>). The above formulas, when expanding the 
exponentials in a loop expansion, possess a representation in terms of Feynman graphs as 
indicated in Fig. [5]. Interestingly, all one particle reducible graphs vanish due to the property 
diG^~° = (c'(0) =0). In addition, since each graph to order U 3 possesses a counterpart 
with a minus sign which is the product of two (factorized) graphs with independent sets of 
loop integrations, this automatically cancels all such (factorized) graphs. 



A. = 

dl 



Figure 3: Graphical representation of the expansion of the modified RG equation in powers 
of the local part up to OiV^ ). The points represent the vertices Ui, the broken solid lines 
are propagators on shell dG, the full solid lines are propagators G, the indices I or fi are 
indicated at each vertex and at each bond. While the solid lines are necessarily present, 
additional dashed lines appear in arbitrary number when performing the loop expansion. 

The RG equation, at any given order in Ui, can be further expanded in the number of 






loops by expanding the exponentials in (|187|) . Let us give the specific result for the case of 



33 



a diagonal gaussian part Gf^ = S^Gf, the generalization being straighforward. To order Uf 
and up to two loops we obtain from (|187|) the RG equation for Ui((j)) as: 





+ + [ d » (^TiM + ifAM + 4 A 'M 

J 


J (189) 


up to 0(Uf) terms, where the contraction graphs are 






Di{4>) = 


dijUiifydijUtic/)) 


(190) 




F,(6) = 


OinmUl ((t>)djimUl(d>) 
^ %j Til «\t J %J ifh w ( \ t / 


(1 Q1 1 




TiM = 


dijUi((j)) djmUp (0) <9 mi ^ (0) 


(192) 




A ( J*\ 

AM) = 


o i jUi(<p)d imn U^)d jmn U^((p) 


(193) 




= 




(194) 


and the momenta 


graphs read 






I? 


= \j q dG 1 G i 




(195) 


If 


1 } A 

A 1 °9l+92 
1 J 9i 9293 




(196) 




= [GfdGfdGl 




(197) 


1 i fJ , 


= / (^91+92 + 
J 919293 




(198) 


I A ' 


J 9l9293 ^ 


(199) 



Note that to two loops the RG flow is generically non local in I. The values of the above 



integrals will be computed in [16 



C.4 RG equation for Ui(<f>) 

The systematic expansion in multilocal interactions can also be performed directly on the 
functional V ((</»). The procedure parallels the previous section and its details are given in 



16]. Here we give only the result for a translationally invariant theory, for the RG flow of 



Ui(4>) in the translationally invariant case to order Uf. It reads: 

SUM) = -^dG^fd^U^) (200) 

1 f r/ X l + X 2\ o12/o^y <ty^y.t=0\ ol2 



JxiX2 £ JO 



4 

1 f Xi + x 2 + x 3 f l 
■7 / <H o ) / d V 



l<fe=</> 



5 123 (9G/ _ dGr )d 123 e -la^ Gl ^ d 12 dG ^3 _ idem Xi = ^ = ^1+^2 

^e-i al2c? - 9l2 9 1 9G^ 2 C/ l ,(0 1 )|^ = ^ l/ (0 2 )|^ = ^ M (03)| 03= ^ 
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Note that the idem x\ 
the G's. We have also defined: 



X2 = Xl + X2 applies only in the square bracket to the arguments of 



d i...p dGd i... P = j2 dG^dfdf 

a,P=l 



(201) 



Again the 1PR diagrams are elimnated by construction since they have one point x« on 
which one can integrate freely, producing a dG q=0 which vanishes by construction. The loop 



expansion of this formula will be detailed in [16 



C.5 Computation of the correlation functions via RG 

The invariance of the generating functional W(J) of the (connected) correlation functions 
with respect to I is now used as a tool for computing the correlation functions of the initial 
model So. The expansion of W(J) in powers of the running interaction Vi(4>) reads formally 
to all orders: 



1 +oo 

W(J) = -J:G:J+J2 K r. 

1 m=l 



(202) 



where we have defined: 



j-.g-.g: 



J (W)) T 



Gi 



Y — ^—^J:G:G l 1 :G:J 



$> m x m = ln(l + r£ 



m=l 



71=1 



-x) 
n\ 



(203) 
(204) 
(205) 



Up to second order, this expansion reduces to: 



W{J) = -J : G : J - e" 



J:G:G7 -G-J 



.,J:G:G 



:0 W) 



i_„-5 J:G-G l :G:J 

+ 2 



e J:G:G-^ Vl{(j)) 2 



(e 2 

g 1 2 V 



Gi 

J-.G-.G^-G-J 



-,J:G:G 



(206) 

^l{<l>)\ G f+0{V!) 



The interaction functional Vi((f>) defined by Q4"2j ) naturally appears in the expansion of 
W(J). Using the properties (J148Q and [«4(0)]g = e2^ :G: ^.4(0), one obtains: 



W(J) = X -J:G: J-V t (G: J) + l(yf(G: J)-\f(G: J))+0(Vf) (207) 

-— ~ _S_.Q ._§_ a 

where Vf(^) = e^i' V;(0i)Vi(02)|^ 1=( /, 2 =^. On this expression, it becomes obvious that 
Vy(J) is indeed /-independent (order by order), as a consequence of the RG equation for 
Vi(4>). As is clear from these formulae, all external legs of correlation functions will carry 
the propagator G while all internal legs will carry Gi. 
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We must now distinguish between the two methods which consist in performing the 
multi-local expansion on V/(0), Vi{4>) respectively. Before doing so, we give a formula, in 
Fourier representation, which is valid in both cases: 

W(J) = \j :G:J - I Ul<e iK - {G - J)x (208) 

2 J xK 

- ( e ^.(G:Jr + iP.(G:J)v (yKP {x _ y) _ l^-KXf? P _ ^(jf) ( 20 9) 
JxyKP V 2 / 

U* = Ufe-^r -* (210) 

V^P( X ) = V K,P( x)e -\K.Gr .K-±P.Gr°.P-K.G*.P (2n) 

the way to compute the functions Ui and Vi being however different in each case. Inserting 
the corresponding formula for Vi as a function of Ui yields expressions in terms of Ui only, 
which we now give in each case (for variety, we also alternate between the -equivalent - field 
and Fourier representations). 

C.5.1 Method with V,(0) 

We start with the formalism using the multi- local expansion of Vi. One finds: 

W(J) = ~J : G : J — I Ui((G : J) x ) + ~ / / [5 IV (e 9 ^~ V 9 * - I) (212) 

2 Jx 2 Jxy Jl'>0 

-6 w d v {e dlG r Vd2 - S(x - y) [ e 9 ^)^,^) | 01=(G:J) ,^(0 2 ) |*,=«w),] (213) 



up to 0(U 3 ) terms. We denote 5w = 5(1 — V) and 9w = 9(1 — /'). From this formula one can 
compute all connected correlations to 0(U 2 ). Let us give the self energy, defined as usual 
from the two point function C = G + 5C as: 

£ = C- 1 - G- 1 = -G-^CG- 1 + G-HCG^dCG- 1 + 0(5C 3 ) (214) 

It reads: 

£«=° = d^fim - I + dld])(e 9lG ' 92 - 1) - d^G**?] U^iPMU^o 

J X 

Yfl - £9=0 = I ^iqx _ ^x ( 215 ) 
Jx 

E?. = -djd] ( i6 u >(e 9lG * 92 - d l G x d 2 - 1) - 9 lv d l dG'td 2 e 9lG ^ 2 }U v {^)U v (h)W^ 
Note that it involves a term with a G propagator. 

C.5.2 Method with V/(</>) 



Inserting the multilocal expansion of Vi, ( |206| ) transforms into an expansion in powers of the 
local interaction Ui(<f)): 

W {J) = - J : G : J - I [ U«e tK - G *-- J + - f f e ^.GMP.G y ,jpKfjP (e -K.Gf-» >.p _ 1} 



2 Jx Jk 2 Jxy Jkp 

f tfUfUfa (e- K - G *> Vp - 5(x -y)f e -K-G?,.P + K.(G t -G°).p\ ] + 0(f/3) (216) 
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Using the RG equation for Ui, it is easily checked again that this expression is /-independent 
order by order. From ( |212| ), one can compute the self energy Tfl of the theory. One gets to 
order Uf: 

= ~ f KiKjUfe-^r * + [ {KiKjAfiq = 0) + K % P 3 B? p {q = 0))(217) 
TP - £«=° = K^Bf^q) - B« p {q = 0)) 

where 

Af p (q) = [ U l «U p e--z KG r° K + PG r 0p (6 l „ l , [ e ^(e- KPG t - 1) (218) 

Jl'>0 Jx 

B« p (q) = I U^U p e^ KG r 0K + PG T 0p {8 l ^ I e ^(e- KG * p + KG X P - 1) (219) 

Jl'>0 Jx 

+6l _ v j ( e <*> _ e K ^- G ^ p )KdG^ v Pe- KG v p ) 



D Dimensional reduction from graphs 
D.l Perturbation theory 

In this appendix, we sketch diagrammatically how the perturbation expansion in R of the 
average of any observable A[u] at T = is the same as the one which would be obtained in 
the Gaussian theory corresponding to a simple random force. 
Precisely, the actions 

<%] = ^= E / - i E / R « - <) (220) 

M a Jxy M ah Jx 

and 

^] = —^—j:[ u* a {{G*G)- l r< (221) 

-Iti (Uj ab Jxy 

where (G * G) x = J Z G G yield the same results when computing the average of any 
functional A[u] of the replicated field at T = 0, e.g. A[u) = T\i u ai- 

To this aim, one first show that the perturbation expansion within ( |220| ) is well-defined 
at T = 0. We use a diagrammatics with propagator 

(u x X) = T5 ab G xy (222) 



which conserves the replica index, and vertex 

-^E j x RK-<) (223) 

associated to one point in space but involving a summation over two replica indices. Thus 
we choose to split the vertex into two subvertices corresponding to each replica index. 
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Figure 4: Typical graph contributing to the computation of the average of an observable A 
to third order (p = 3) in perturbation theory. Note here that v = 2p — 6, t = 1, k = 3, 
e = 4, c = 3, / = 1. 



For any graph occuring in the computation of (A[u]R p )g, let us denote by k the number of 
lines connecting A to the vertices (involving the extraction of k legs from A : d u *i M *kA[u\). 
Let /C be the graph obtained by considering only the splitted vertices, the propagators 
between them, forgetting the observable and the k lines attached to it. The graph K, has v = 
2p subvertices. Contrarily to the initial graph with unsplitted vertices, K, is not necessarily 
connected and is made of c connected components. To each one corresponds a replica index. 
If one of them is not connected to the observable, i.e. if it does not inherit from a replica 
index contained in A, then the summation over this index is free, giving a factor n. Hence 
each connected component has to be linked to the observable in order to survive the n — > 
limit, which yields k > c. 

Collecting the factors of T in front of the initial graph (/C, A, the k propagators, and 
possibly t tadpoles on A), the power ofT is t + e + k — v where e is the number of propagators 
in /C. Euler relation in K, reads v + I = c + e where I is the number of loops in /C. But since 
k > c and I > 0,t > 0, one obtains that each graph is in factor of a non-negative power of 
T. The existence of the T = perturbation theory is thus confirmed. 




Figure 5: Vanishing contribution to the average of A : the second graph is obtained by 
mounting a branch to the first graph. The only remaining graphs are simple "flowers" . 

The graphs which remain at T = have t — 0, I — 0, k — c, which means the following 
properties (i) their subgraph /C has no loop, each of its connected component is a tree, (ii) 
there is no tadpole on A, (iii) A is linked to each connected component of /C by one unique 
propagator. This result is easily extended to a non-Gaussian disorder, which possess higher 
cumulants of the general form (pU|). 

The second part of the argument uses the property of translation invariance in u space 
of the disorder distribution, on which the first part does not rely. Since each connected 
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component of /C is linked to A, let us call root the point to which it is attached. This 
provides a natural orientation to the branches of the trees from root to leaves. If any point 
of /C possess at least a branch going to the direction of the leaves, the graph obtained 
by mounting this branch to the companion point (which belongs to the same R before 
splitting) has the opposite value. One can convince oneself that such graphs can be grouped 
by mutually cancelling pairs. Thus the only graphs which survive to this mounting operation 
look like flowers, with A at the center and petals R made of two propagators (see Fig. [5]). 

The generalization including higher cumulants is straightforward, but yields a non- 
Gaussian theory. The corresponding equivalent action for computing observables is 

s r fM = ^E/ <(g">^-E^^...^ (7V) (o...o)/(E<) 7V (224) 

Al a Jx V N>2 iV !i Jx a 

Even if this is not obvious on (p24|) , this action possess statistical tilt symmetry as can be 
checked thanks to n — > 0. 



D.2 Corrections to R 

The computation of the effective action T (1PI) involves corrections to the various cumulants 
of the disorder. At T = 0, the graphs correcting a N replica term (iV th cumulant) is made of 
N connected components, so that there exists a free sum over N replica indices. The power 
of T in front of such a cumulant has to be — N. The graphs correcting R with p i?'s are 
made of two connected components (c = 2), the power of T is e — v where e is the number 
of propagators and v = 2p. Euler relation yields e — v = I — c > —2 with equality for I = 0. 
Hence such graphs are made of two trees. 

Furthermore, the graphs such that the two points of a splitted R are connected to the 
same connected component and such that one of them is connected to at least two branches 
vanish. This can be seen by mounting one of these two branches on the companion point. 
Hence, if two points of a R belong to the same connected component, then each one is 
connected to a unique branch. As a corrolary, the two points of an R cannot be connected 
to each other by a branch (since it would be impossible to connect this R to the rest of the 
diagram thanks to the argument above). 

This considerably reduces the form of the possible corrections to R. These corrections 
obey in particular 

5R=(e-40R + (uR' + ^—j RP (225) 

where the last term symbolically only means that the p th order term contains 4(p — 1) 
derivatives (not that it is a total 4(p — 1) derivative). We allowed for a field rescaling with 
exponent (. To order R 3 the arguments above allow only for the following corrections 

5R = (e- 4()R + CuR' + K(^R" 2 - R"R"(0)) + A(R" - R"(0))R"' 2 + C(R" - R" '(O)) 2 #"(226) 

with some constants K, A and C, and valid only for an analytic R(u). In the periodic case 
(C = 0), the fixed point equation is easily solved since there exists to any order in e a fixed 
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point function of the form 

R*(u) = a + bu(l -u) +c{u{l -u)f (227) 

where a, b, c can be computed in series of e, once the coefficients of the fixed point equation 
are known. This is further examined in Sections IR 



E variational calculation 

Here we sketch the derivation of the scaling function for the confined interface using the 



replica variational method, extending the explicit solution of Ref. [27] to a non zero mass. 
We use all notations of Ref. [27] and Ref. |23| . 



The disorder correlator for the random field problem studied here corresponds to the 
case 7 = 1/2 and g = a for the parameters of E7J. Applying the variational ansatz for 



N = 1 components yields the function f(x) = gJr 2 + x which describes the correlations 



and g = ay 2/ it. We have artifially extended the correlator to small scales, so as to obtain a 
well defined T = limit. The large scale results however are independent of the small scale 
details in the limit of small m. The variational equations reads: 

(rj + 2TB(u)f/ 2 = gj D ([a)(u) + m 2 )^ 2 (228) 

with jjj = J k (k 2 + 1)~ 2 = T(2 — d/2) / ' (Ait yd/2) and the equation for the breakpoint (see Ref. 
PI) is: (rj + 2TB(u c )f /2 = gj D (Z c + m 2 )-^ 2 . This yields the solution: 

[a] (u) +m = Au 2/e = Z c (u/u c ) 2/0 u* <u<u c (229) 

a(u) = a(0) = A e / 2 -^—m 2 - 6 u < u* (230) 
2 — 

and [#"](«) = S c for u > u c . Here 6 = is the free energy fluctuation exponent and A = 

(-§fg 2 ^j]j^) 2 ^- This solution allows to compute the small q behaviour of the correlation 
(for small m ~ q) asL 

T r-c du K(u/u c ) 2 / d - m 2 1 

U « U -« q 2 + m 2 ^ + 7 Uc(m2/Ec)e/2 U 2 q 2 + T, c (u/u c fl e K ' q 2 + m 2) 1 ' 



which yields the large scale result 143 given in the text. In D = one recovers the result of 
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u 2 = ^m-^ (232) 



F FRG to two loops 
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F.l Method with V 



The exact RG equation to order U 3 given in |C4| when expanded to two loops yields the 
following finite temperature RG equation for Ri(u) at large I (the derivation is detailed in 



m = eR l + ffi; + K lfl ^{uf - jjj(0) JiJ(u)J (233) 

+K A (R" - R"(0))R"' 2 + K C \R" - R"(0)) 2 R"" (234) 
+K s d 12 S l (u, u, 0) + T^l (R'?(R;(u) - ^(0)) - R'^m;) + T^; 2 (235) 

In this formula all terms of order R 2 and higher are retarded, and integrals J d/i are un- 
derstood. For the R 3 terms we have ommitted the retardation integrals (which involve an 
additional integral ftf dp) because near the fixed point they can be replaced by a single num- 
ber. The feedback of the three replica term is through its partial derivatives. This three 
replica term satisfies its own RG equation given in |16|]. The precise values of all coefficients 



are detailed in |[T^ 



To obtain the large I limit of this equation and thus the fixed point equation we use the 
fact that Ti — > and the property (|1 1 7|) which we need only to lowest order. Integrating out 
the third cumulant RG equation (which reaches the fixed point value ( |103| )) and inserting 
into (|235|) yields the FRG fixed point equation (for simplicity in the periodic case): 

= e R* + ^R*" 2 - R*"(0)R*" (236) 

+ (- + K°)(R*" - R*"(0))R*"' 2 - -R*"'(0 + ) 2 R*" + K° (R*" - R*"(0)) 2 R*"" (237) 
2 2 

(j 

where only one coefficient, K , is non universal. We have absorbed the exact one loop 
coefficient K obtained in (|9l"l) in R (whose first 0(e) correction is non universal). In terms 
of the cutoff function c(s) = J a c(a)e~ as (or its expression as a sum of exponentials) one has: 

K = 2(2!) rf / 2 [1 + ^2 ~ la £(G) llla + lL £(fl)£(6) lll(a + 

K° = - [ 6(a) In a + / c(a)c(6) ln(a + b) = f° ds ~ c< ^ s >> (238) 

J a Jab JO S 

to the desired order in e. This equation has a fixed point solution: 

R*(0) = K^H^ " ^« 2 (l - nf) + (239) 

2 A + 10K° K° n , 1 + 3K° 2/1 . 2 ,. /0 . n , 

e( — === u(l-u) u 2 (l-uY)} (240) 

v 7776 216 v ; 108 v ; n v ' 

Note that this solution, for K =0 (i.e hard cutoff) does not contain a term with R'(0 + ) ^ 
( "supercusp" ) 
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F.2 method with V 

The method using the Wick ordered functional can also be used. Using the RG equation 



to third order in U and two loops ( |189| ) one finds the equation for the corresponding two 
replica part of U : 

m = (e - 4C)i* + CuR't + K Qi?,» 2 - R!{ '(0) R? '(«)) (241) 
+K A (R" - R"(0))R m + TiKfR"' 2 (242) 



with: 



K = ~IFM (243) 



where the integrals have been defined in Section |C.3| . Calculation shows that the constant 
K A = 1/2 (to lowest order in e) independent of the cutoff function c(s). The coefficient of 
this term is in agreement with [0|. The analysis of the boundary layer is more intricate 
in this formulation 
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